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A method of targeting engineering synchronization states in two identical and mismatch 
chaotic systems is explained in details. The method is proposed using linear feedback 
controller coupling for engineering synchronization such as mixed synchronization, linear 
and nonlinear generalized synchronization and targeting fixed point. The general form of 
coupling design to target any desire synchronization state under unidirectional coupling 
with the help of Lyapunov function stability theory is derived analytically. A scaling 
factor is introduced in the coupling definition to smooth control without any loss of 
synchrony. Numerical results are done on two mismatch Lorenz systems and two identical 
Sprott oscillators. 
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1. Introduction 

Chaotic system is defined by its complex dynamical behaviors, particularly their 
extreme sensibility on initial conditions and parameter variations, which make their 
behaviors long-term unpredictable. The coupled chaotic system exhibit interesting 
complex feature like riddled basin, on-off intermittency, different types of synchro¬ 
nization and amplitude death between the oscillators. In 1665, Huygens discovered 
an odd kind of sympathy in two pendulum clocks suspended on a beam. This was 
visualized as an interesting example of synchronization phenomena between two 
clocks. Recently in last two decades this phenomenon has become an important 
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topic to be discussed in the aspect of nonlinear coupled systems i.e. mainly in 
chaotic system. 

Pecora and Carroff introduced a method to synchronize two identi¬ 
cal systems with different initial conditions. Several types of synchroniza¬ 
tion like complete synchronization(Csj21, antisynchronization ( ASj^l^lSI, lag 
synchronizationt'LsI-' 7 , phase synch ronization (PSjpl, generalized 

synchronization(GSprojective syn chr oniz ation an alternate ap¬ 

proach like engineering synchronization^ ^ II 2Q] e ^ c have been reported before. 
In CS the state variables of coupled system are completely coincide where as phase 
difference 7r exist between them for AS. The state variables of the coupled oscil¬ 
lators emerge into a mixed synchronization (MSwhere a pair of state 
variables develops a CS state while another pair is in a AS state. It is true that exis¬ 
tence of a MS state was reported earlier^ in coupled co-rotating Lorenz systems for 
a specific scalar coupling. This emergent MS sta te has justification in the inherent 
axial symmetry of the chaotic flow and majP21. Recently MS sta te is emerges in 
counter- rotati ng oscillators under linear diffusive couplin^SElES] wag recently 
reported^® 23 that, given a definition of a model system, one can design a coupling 
to target MS state in delayed and nondelayed systems. The combination of CS and 
AS provide a way to cryptographic encoding for digital signal through parameter 

modulatioif^ 2S1 

The other kind of synchronization namely GS was introduced which tried to ex¬ 
plain emergence of a kind of functional relationship between the coupled oscillators 
in drive - response mode. GS is thus seen as an evolution of a functional relation¬ 
ship between state variables of driver-response system. This functional relationship 
for GS is generally unknown under conventional diffusive coupling. Ano ther most 
intriguing effects in dynamical system is so-called amplitude death (AD)^ 3H23 
dealing with suppression of oscillations to a steady state. Amplitude death in cou¬ 
pled systems generally occurs when the amplitude of the oscillation is damped out to 
a steady state. This all types of dynamical phenomenon also reported earlier using 
different coupling for getting desired state. For getting desired synchronizatio n state 
different coupling configuration has been investigated such as unidirectionaPSESl 
bidirectional^, repulsive^®, inhibitory^ ^8 or excitatorjEH, synaptic^coupling etc. 
To observe different types of synchronization state between two or more oscillators, 
major efforts are given on the strength of coupling between them and parame¬ 
ter mismatch or noise intensity. Under linear coupling, in most studies, different 
synchronization state i.e. CS, AS, PS, LS, GS are observed by varying the cou¬ 
pling strength but above critical coupling strength most of cases desynchronization 
regime observed. 

Recently the concept of engineering synchronization^ in nonlinear oscillators 
given importance, for practical application. In engineering synchronization, it is as¬ 
sumed that the definition of chaotic flow is known but the definition of coupling 
unknown apriori. In this paper we address the question “how to define the cou¬ 
pling between two coupled oscillators to get desired synchronization state when the 
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definition of flow is known?”. 

In this paper, we propose linear feedback controller (LFC) to target engineering 
synchronization state. We design coupling definitions to target engineering mixed 
synchronization, linear and nonlinear generalized synchronization and target fixed 
point. Partial symmetry of the system is the necessary condition for MS state. Note 
that this restriction is valid for MS under linear conventional diffusive coupling too. 
Here, we remove this restriction to get MS state. Given a model system and target 
synchronization state, the coupling is first defined analytically to establish engineer¬ 
ing synchronization state between two chaotic systems. The stability condition for 
targeting engineering synchronization state is derived using Lyapunov functional 
theory. A smooth control from one synchronization state to another synchroniza¬ 
tion state namely CS to AS and vice-versa is achieved by inserting a scaling matrix 
in the definition of coupling. This is the main contrast to what usually observed 
for conventional linear diffusive coupling. Engineering linear and nonlinear gener¬ 
alized synchronization is also observed for a suitable choice of target function. In 
GS state the functional relationship between drive-response systems is generally 
unknown under conventional diffusive coupling. But targeting engineering general¬ 
ize synchronization, the functional relationship is known which is a contrast to the 
conventional diffusive coupling. We also target the fixed point which may be a fixed 
point of the system or any other fixed point whereas the driving system is in chaotic 
state and response system target to a fixed point. Amplification and attenuation of 
response system’s attractor is also possible using LFC coupling. The proposed LFC 
coupling is not dependent on system’s parameters. For a given system model and 
a desired synchronization state, if the coupling term designed then smooth control 
of one synchronization state to another synchronization state with amplification or 
attenuation is possible without loss of synchrony. The scaling of a chaotic attractor 
has potential application in secure communication. Numerical simulations are done 
using mismatch Lorenz system and identical Sprott oscillators. 

We organize the paper as follows: Sec. II, general theory of linear feedback con¬ 
troller (LFC) is discussed for targeting engineering synchronization between two 
non-identical systems. The theory is also valid for identical systems taking mis¬ 
match term as zero. In Sec.Ill, the theory for LFC using unidirectional coupling is 
elaborated how to realize MS state i.e. co-existence of CS and AS between differ¬ 
ent state variables using numerical example of non-identical Lorenz system®. Nu¬ 
merical example on identical Sprott oscillator — for engineering linear generalized 
synchronization is discussed in Sec. IVA. Sec. IVB is elaborated with engineering 
nonlinear generalized synchronization using identical Sprott oscillator. In Sec. V, 
targeting fixed point is discussed on Sprott system where the driving states are in 
chaotic state and response states are in a desired fixed point state. The results are 
summarized in Sec. VI. 
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2. General theory of coupling 

We consider the driver system as 

x = f{x,n) + Af(x,n) (1) 

where x £ R n , y is the vector of parameters and A f(x, y) = f(x, /i + 6y) — f(x, y), 
in general and Sy denotes the mismatch parameters. Otherwise if all parameters 
are appeared in linear form in /(•) then A f(x,y) = f(x, Sy). 

Now we consider the response system as follows 

y = f(y,v) + u (2) 

where U is controller which we want to derive using Lyapunov stability theory. The 
error signal of the coupled system is defined as 

e = y-</>(x) (3) 

where is the target or goal function. Error dynamics is 

e = y- D<j>(x)x = f(y, /x) - D(j>{x)[f{x, /x) + A f(x, /x)] + U (4) 

where D<f>(x) is the Jacobian matrix of the function <p(x). We choose the controller 

U as follows 


U(t) = u'(t) + w(t) (5) 

where u'(t) and w(t) are the active and linear feedback controller respectively. So the 
error system to be controlled is now a linear system with the control input function 
w(t) as function of the error states. When the error system will be stabilized by the 
feedback w(t), the error will converge to zero as t —> oo which implies that the drive 
and response system are globally synchronized. To achieve this goal we choose w(t) 
such that, 


w(t) T = A.e(t ) 

where A is a matrix of order n x n. So A should be chosen properly to get target 
synchronization state. Now we consider the Lyapunov function as follows 

V = ie T e (6) 

where T is transposition of matrix. It can be easily verified that V(t) is a non¬ 
negative function. Assuming the system’s parameters are all known, we make an 
appropriate choice of the controllers u'(t) and w(t) so that ^- < 0. The error 
dynamics will be asymptotically globally stable if ^- < 0 and thereby realize any 
targeted synchronization state. 
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3. Targeting Engineering Mixed Synchronization 

Consider the three dimensional autonomous chaotic Lorenz system^, described by 
the set of equations: 


X\ = <r(x2 — X\) 

X2 = Tx 1 - X 2 — X\Xz 

x 3 = X 1 X 2 - fix 3 (7) 

where xi,x 2 ,x 3 are the state variables and a,r,/3 are parameters of the above 
system. The system exhibits chaos when a = 10, /3 = |, r = 28. 

We consider the system (7) as the driving system with parameter mismatch as 
the following equation ( 8 ) where 5a, 5r, 5f3 are the parameter mismatch correspond¬ 
ing to a, r, /? respectively. 


Xi = {a + 5 a){x 2 - x\) 
x 2 = (r + Sr) x\ — x 2 — X1X3 

±3 = X\X 2 - (P + 5(3)x 3 ( 8 ) 

We consider the response system as follows 

Vi = 0 -(y 2 ~Vi) + »i 

V2 = ryi - V2 ~ 2 / i 2/3 + u 2 

2/3 = 2 /i ?/2 - PV3 + u 3 (9) 

where u\, u 2 , u 3 are the controllers to be chosen later. Here our aim is to 
determine the controller for the purpose of mixed synchronization with parameter 
mismatch. Let the target state is defined as 


tp(x) = ax (10) 

where a = ( aij) n xn is the scaling matrix. For mixed synchronization we choose the 
scaling matrix as a = diag(a 11 , a 2 2 i 0 : 33 )- When an = 1 ,i = 1,2,3, we can only 
achieve CS state without any amplification or attenuation. For \au\ > 1 or \aa\ < 
1, i = 1,2,3, scaling of the size (amplification or attenuation) of the driving attractor 
at the response attractor is also possible. 

The error system is defined by 

ei = 2/1 — 0111X1 
£2 = 2/2 — 0:222:2 

e 3 = 2/3 - 033 X 3 ( 11 ) 

where 011 , 022,033 are the scaling factors. Then the error dynamics can be 
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written as 

ei = cr(y 2 - yi) + «i - an(cr + Aa)(x 2 - Xi) 

= er(e 2 — ei) + cra 22 x 2 + Ui — au<rx 2 — an A a(x 2 — Xi) 

= a(e 2 - ei) + ax 2 (a 22 — an) — anAa(x 2 — x\) + u\ 

= a(e 2 - ei) + w\ ( 12 ) 

where u\ = an Aa(x 2 — x\) — a (a 22 — an)i2 + Wi 

<32 = ryi - 2/2 - 2/12/3 + u 2 - a 22 (r + Ar)xi + a 22 x 2 + 
ot 22 x\x 3 

= r(y\ - anii) + rauii -e 2 - 1/12/3 + u 2 
—a 2 2 [r + Ar)x\ + a 22 x 3 x 3 

= rei — e 2 + w 2 (13) 

where u 2 = 1/12/3 — a 2 2XiX3 + a 2 2(r + Ar)ii — rauii + w 2 


e 3 = 2 / 12/2 - Py 3 + m 3 - 033 X 10:2 + a 33 (/3 + A/3)x 3 

= —/?e 3 + W 3 (14) 

where u 3 = - 2 / 12/2 + 033 X 1 X 2 - a 33 A/3x 3 + u > 3 

Thus, the system (12 — 14) to be controlled is a linear system with the control 
input function w(t) = [wi(t),w 2 (t),w 3 (t)] T as functions of the error states. When 
system (12 — 14) is stabilized by the feedback w(t), the error will converge to zero 
as t — > 00 which implies that the system (8) and (9) are globally synchronized. To 
achieve this goal, we choose w(t) such that, 

[w 1 (t),w 2 {t),w 3 (t)] T = A[ei{t),e 2 (t),e 3 (t)] T (15) 

where A = (a,y)3x3 is any matrix. Different choices of the matrix A are possible. 
We make a choice of the matrix A is as follows: 

/ 0 -cr 0\ 

A= -r 0 0 (16) 

V 0 0 0/ 

Other choices of the matrix A are always possible. We re-define the feedback control 
function as follows: 


/ Wi \ / 0 — erO\ / e \ \ 

\ w 2 = —r 0 0 I x e 2 

\ w 3 ) \ 0 0 0/ \ e 3 / 


( 17 ) 
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Fig. 1. For mismatch Lorenz system ( a ) time series of x\ in blue line and atunuated y\ in red 
line with scaling factor an = 0.5 in CS state, (b) represents time series of X 2 in blue line and y 2 
in red line with c *22 = — 1 in AS state, (c) time series of xs in blue line and response ys in red 
line in CS when Q 33 = 1. and (d), (e), (f) are corresponding synchronization plot of (a), (b), (c) 
respectively. The parameters mismatch are A a = 1.0, A r = 0.0, A/3 = 0.0. 


Hence the error system becomes 

e’i = — crei 

e ‘2 = -e2 

e 3 = -/?e 3 (18) 

Now we consider the Lyapunov function as follows: 

V (t) =-(el + el + e l) (19) 

By using control law, the time derivative of V is obtained as 

V ( t ) = eie'i + e 2 e 2 + e 3 e 3 

= ei(-aei) + e 2 (-e 2 ) + e 3 (-/3e 3 ) 

= -ae\ -4- pe 2 3 (20) 

Now from the Lyapunov stability function, we can say that the error system 
equation is asymptotically stable if V < 0, if cr > 0 and /? > 0. We integrated 
numerically system (8) and (9) with the controllers (iti,u 2 ,u 3 ) using fifth order 
Runge-Kutta-Fehlberg method with integration step At = 0.001 and taking ran¬ 
dom initial conditions. Numerical results of MS state are shown in Fig. 1 for two 
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coupled mismatch Lorenz oscillators eqn. ( 8 ) and eqn. (9). The time series for three 
pair of variables (aq, yi),(x2, 2/2) and ( 2 : 3 , y 3 ) are plotted in Figs. 1(a), 1(b) and 
1(c) respectively. Synchronization plot also shown in Figs. 1(d), 1(e) and 1(f) where 
(xi,y 1 ) and (£ 3 , 2 / 3 ) are in CS state whereas ( 2 : 2 , 2 / 2 ) is in AS state. This MS state 
with amplification or attenuation is not possible under linear conventional diffu¬ 
sive coupling between two non-identical chaotic systems'^. This is contrast of LFC 
coupling with the diffusive coupling. 

4. Targeting Generalized Synchronization 

A generalized synchronization scenario was reported earlier using linear diffusive 
coupling between two nonidentical or large parameter mismatch oscillators. In GS 
scenario the functional relationship is usually unknown. To study line ar g eneralized 
synchronization, we consider the three dimensional Sprott oscillator^!, described 
by the following sets of equations: 

*1 = X\X2 — *3 
*2 = Xi — X 2 

*3 = Xi + ax 3 (21) 

where X\,X2,x 3 are the state variables and a is the only parameter of the above 
system. The system exhibits chaos when a = 0.3. 

The system (21) with the controllers u 1 , M 2 , M 3 to be chosen as driving system 

*1 = X\X2 — x 3 + Mi 
*2 = *1 - x 2 + m 2 

*3 = *1 + ax 3 + M3 (22) 

So the response system as follows 

2/1 = 2/12/2 - 2/3 
2/2 = 2/1 - 2/2 

2/3 = 2/i + a//3 (23) 


4.1. Targeting Linear Generalized Synchronization 

For linear generalized synchronization, we consider the target function as 


x) = 



= ax = 


an OU 2 a 13 
a 2 l Q?22 Qf23 
,0:31 0:32 0:33 


\ , 

f xi ) 

r 

* 2 

/ 1 

^*3/ 


(24) 


Now our aim is to determine the controller for the purpose of linear generalized 
synchronization. Let the error vectors are defined as 


ei= 2 /i - Oil 1X1 — ai2X2 — ai 3 x 3 
02 = 2/2 — 0:21*1 — 022*2 — 023*3 
£3=2/3“ O31CC1 — 032*2 “ 0332:3 


( 25 ) 
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where (o^)3 X 3, ( i,j = 1 , 2 , 3 ) are the scaling matrix. 

Then error dynamics can be written as, 

ei = 2/12/2 - 2/3 + «i - an ( xix 2 - x 3 ) 

-<212(2:1 - x 2 ) - <213(2:1 + <<2:3) 

= 2/12/2 = ui - anxix 2 - ( y 3 - <2312:1 
—<2322:2 — <2332:3) — <2310:1 — 0322:2 
-<2332:3 + <2112:3 - 012(2:1 - x 2 ) - 013(2:1 + ax 3 ) 

= - e 3 + wi ( 26 ) 

where «i = 0112:12:2 - 2/12/2 + 2:1(031 +012 + o i3 ) + 2:2(032 - oi 2 ) +2:3(033 - o n + 
0013) + Wi 

e 2 = 2/i - 2/2 + u 2 - 021(2:12:2 - 2:3) - 022(2:1 - x 2 ) 

-023(2:1 - 0x3) 

= ei — e2 + on2;i + 0122:2 + 0132:3 — 0212:1 
—0222:2 — o 232:3 + u 2 — 0212:12:2 + 0212:3 — 0222:1 
+0222:2 — 0232:1 — 00232:3 
= ei — e 2 + w 2 

where u 2 = ii(-Qn + o 2 i + a 2 2 + 023) + X2(—a 2 i) + x 3 (—ai 3 + a 2 3 
0023) + 0212:12:2 + w 2 

e 3 = yi + ay 3 + u 3 - 031(2:12:2 - 2:3) - 032(2:1 - x 2 ) 

-033(2:1 + 02:3) 

= ei + ae 3 + (onXi + 0122:2 + 0132:3) + 0(0312:1 
+0322:2) + 0332:3 + u 3 — 031X32:2 + 

031X3 - a 32 (xi - x 2 ) - a 33 (xi + ax 3 ) 

— s ± + ae3 + w 3 ( 28 ) 

where u 3 = xi(—an — aa 3 i + 031 + 033) + X2(—012 — 0032 — 032) + X3 — 2031 + 
a 3 ixix 2 +u> 3 

Particular choice of linear feedback control function as, 

ei 
e2 

e 3 

Hence the error system become, 

ei = -ei 

e 2 = —e 2 



/-I 0 1 
= 1—10 0 
\-l 0 —2a 


(27) 
— 021 + 


e 3 = -ae 3 


(30) 
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Fig. 2. For Sprott system (a) time series of x\ in blue line and y\ in red line, (b) represents 
plot of x\ with yi. (c) synchronization mainfold of (yi,0i), (d) (^ 2 ,^ 2 ) and (e) (ys, ^ 3 ) where 
4>l = anxi + a 12^2 + &13 x 3i <f> 2 = <*21 + Ot-22 x 2 + &23 x 3y 4>3 = «31^1 + &32 x 2 + <^33 x 3- All 

scaling parameters value are Qjj = 1, i,j = 1, 2,3. 


The Lyapunov stability function 



-{e\ +4 + ae\) 


(31) 


System is stable because V < 0. The time series, phase portrait, and synchro¬ 
nization plot of identical Sprott system is shown in Fig. 2. In Fig. 2(a) shows time 
series of x\ (blue color) and y\ (red color) and which apparently shows no correla¬ 
tion by visual check. The synchronization plot in {x\, y\) plane appeared to show no 
correlation. However, y\ vs. <fi i in Fig. 2(c) shows clearly 1:1 correlation and confirm 
that the targeting linear GS state. Other linear GS plot in ( 2 / 2 ^ 2 ) and ( 2 / 3 , <(> 3 ) are 
shown in Figs. 2(d) and 2(e) respectively. 
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4.2. Targeting Nonlinear Generalized Synchronization 

For nonlinear generalized synchronization we have taken a nonlinear target function 
such as, 

( x \\ 

<p{x)=\xl\ (32) 

The synchronization errors become 

ei = 2 /i — x\ 
e 2 =U 2 - x\ 

e 3 =y 3 - x 3 (33) 

So the error dynamics as follow 

ei = j/i — 2xi(±i) 

= 2/12/2 - 2/3 + mi - 2xi(xix 2 - x 3 ) 

= 2 / 12/2 + Ml - 2 ^ia ;2 + 2x 3 x 3 - (y 3 - x 3 ) - a : 3 
= —e 3 + W! (34) 

where u 3 = — 2 / 12/2 + 2 a;i(a;ia ;2 — a; 3 ) + x 3 + w 3 


e-2 = V2- 2 x 2 (x 2 ) 

= 2/1 - 2/2 + m 2 - 2x 2 (a; 1 - x 2 ) 

= ( 2/1 - ®i) + x\ - ( 1/2 - * 2 ) + x\ - 2x\x 2 + m 2 
= ei — e 2 + w 2 (35) 


where u 2 = 2 x 1 X 2 — — xf + w 2 


e 3 =y 3 ~ x 3 

= 2/1 + «2/3 + m 3 + xi - ax 3 
= ( 2/1 - x\) + a(y 3 - x 3 ) + u 3 + x\ - x\ 

— e\ + ae 3 + w 3 ( 36 ) 

where u 3 = —x\ + x\ + w 3 

The feedback controller for nonlinear generalized synchronization, we choose as 


( W\ 

w 2 

W 3 


/- 10 1 \ 
-10 0 x 
V 1 0-2 a) 



( 37 ) 
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y. 




y 2 y 3 


Fig. 3. For Sprott system nonlinear generalized synchronization plot of (a) (yi,xi), (b) ( 3 / 2 ,£ 2 ) 
and (c) ( 133 , 3 : 3 ) where an = 022 = a 33 = 1 . 0 . 


Hence the error system become, 


ei = —ei 


e 2 = -e 2 


e3 = —ae3 


(38) 


The Lyapunov stability function 


V — ~( e l + e l + eg) — — (ef + ej + aeg) < 0 


(39) 


Fig. 3 shows nonlinear generalized synchronization plot. Figs. 3(a), 3(b) and 3(c) 
shown synchronization plot in (y 1, Xi), (y 2, x 2 ) and (2/3, x 3 ) respectively. By checking 
it is confirmed the onset of GS and satisfying the relations y 1 = xf, y 2 = x\ and 
2/3 = x 3 . 
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5. Targeting fixed point 

For targeting fixed point of response system the earlier Sprott system(21) is con¬ 
sidered and target function as 

/ a n 0 0 \ 

ip(x) = a = I 0 o 2 2 0 J (40) 

\ 0 0 033 / 

The error sytem comes 

ei = 2/i — an 
e 2 = 2/2 ~ « 22 

e3 = 2/3 — «33 (41) 

So the error dynamics as follows 

ei = jfi = 2/U/2 - 2/3 + mi 

= 2/12/2 - (2/3 ^ a33) ^ 2/3«33 + Ml 
= -e 3 + Mil (42) 

where Ui = 0331/3 - 1/12/2 + wi 

e 2 = 2/2 = 2/i - 2/2 + m 2 

= (2/1 — on) + u 2 + on — o 22 

= ei - e 2 + ic 2 (43) 

where w 2 = —on + o 22 + w 2 


e3 = Vz = 2/1 + ^2/2 + m 3 

= (2/1 - on) + a(2/3 - o 33 ) + u 3 + on + o«2 2 
= ei + ae3 + W 3 


where 113 = —an + 0033. We choose the linear feedback function as 


( W\ 

w 2 

M’3 


/-!0 0 \ /d\ 

-10 0 x e 2 
V 0 0 -2a/ \ e 3 / 


Finally the error system becomes 


ei = —e\ - e 3 


e 2 = —e 2 


d = ei - ae 3 


(44) 


(45) 


(46) 


So ei = e 2 = e 3 = 0 is asymptotically stable using Lyapunov stability theory. 

It is also interesting that AD can be induced in the response system by taking 
on = —1.0, o 22 = 0 . 1,033 = 0.5 as shown in Fig. 4. The driver x(t) is oscil¬ 
latory while the response y(t) ceases to oscillate and is stable at targeted state. 
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Fig. 4. Time series of (a) (# 1 , 2 / 1 ), (b) (# 2 , 2 / 2 ) and (c) (# 3 , 2 / 3 ) shows the driving system states 
(# 1 , # 2 , # 3 ) (blue color) are in chaotic state and response system states ( 2 / 1 , 2 / 2 ,2/3) (red color) are 
in targeted fixed point state. Here an = —1.0, a 22 = 0 . 1,(233 = 0.5. 


The coupling actually stabilizes the response system to any desired state. For 
an = 0-22 = <233 = 0, the response system stabilized at origin (figures not shown). If 
the response system has no equilibrium point at origin, the coupling has an inherent 
property to create a new equilibrium at origin and stabilize the response system at 
origin. For an = —1.0, 0:22 = 0.1, and 033 = 0.5, the response system stabilize at 
yi = —1.0, 2/2 = 0.1, 2/3 = 0.5 are shown in Fig. 4(a), 4(b), and 4(c) respectively 
plotted by solid red color lines. Whereas the driver states xi(t),X 2 (t) and X 3 (t) are 
in chaotic state in Fig. 4(a), 4(b), 4(c) plotted by solid blue color lines. This is one 
of the major advantages of this coupling over the conventional diffusive coupling. 
This targeting fixed point is used in practical system where chaos is an unwanted 
event. 


6. Conclusion 

We explore linear feedback controller based coupling design using Lyapunov func¬ 
tion stability theory for targeting engineering synchronization in chaotic oscillators. 
The general scheme for coupling design in an unidirectional couple oscillators is 
discussed. By introducing a scaling factor in the definition of coupling that al¬ 
lows amplification or attenuation of one attractor relative to another. We described 
the theoretical details of the method about how to design coupling and illustrated 
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with numerical examples o f Lorenz and Sprott systems. Compared to previous cou- 
pling — 34 35,36 37 38|39| 4(3 our ji near feedback controller based coupling has the 
following advantages: (i) previously mixed synchronization is observed using diffu¬ 
sive coupling in a chaotic flow with partial symmetry. We remove this restriction 
for mixed synchronization, (ii) Smooth control from one synchronization state to 
another synchronization state is possible by changing the scaling matrix in the cou¬ 
pling definition, (iii) Linear feedback controller coupling is independent on system 
parameters, (iv) In targetting engineering generalized synchronization state, the 
functional relation between drive and response state is known. But in conventional 
diffusive coupling, this functional relation for generalized synchronization in un¬ 
known always, (v) Amplification or attenuation of response system’s attractors is 
possible. The scaling of chaotic attractor increases sucurity in chaos cryptography, 
(vi) Targeting engineering linear and nonlinear generalized synchronization using 
linear feedback coupling satisfying Lyapunov function stability theory were not ex¬ 
plored earlier. Physical realization of engineering synchronization using electronic 
circuit is also a challenging problem, will be the future work. 
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